Functional renormalization group for anisotropic depinning and relation to branching 

processes 
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Using the functional renormalization group, we study the depinning of elastic objects in presence of 
anisotropy. We explicitly demonstrate how the KPZ-term is always generated, even in the limit of vanishing 
velocity, except where excluded by symmetry. This mechanism has two steps: First a non-analytic disorder- 
distribution is generated under renormalization beyond the Larkin-length. This non-analyticity then generates 
the KPZ-term. We compute the /3-function to one loop taking properly into account the non-analyticity. This 
gives rise to additional terms, missed in earlier studies. A crucial question is whether the non-renormalization 
of the KPZ-coupling found at 1-loop order extends beyond the leading one. Using a Cole-Hopf-transformed 
theory we argue that it is indeed uncorrected to all orders. The resulting flow-equations describe a variety of 
physical situations: We study manifolds in periodic disorder, relevant for charge density waves, as well as in 
non-periodic disorder. Further the elasticity of the manifold can either be short-range (SR) or long-range (LR). 
A careful analysis of the flow yields several non-trivial fixed points. All these fixed points are transient since 
they possess one unstable direction towards a runaway flow, which leaves open the question of the upper critical 
dimension. The runaway flow is dominated by a Landau-ghost-mode. For LR elasticity, relevant for contact 
line depinning, we show that there are two phases depending on the strength of the KPZ coupling. For SR 
elasticity, using the Cole-Hopf transformed theory we identify a non-trivial 3-dimensional subspace which is 
invariant to all orders and contains all above fixed points as well as the Landau-mode. It belongs to a class 
of theories which describe branching and reaction-diffusion processes, of which some have been mapped onto 
directed percolation. 



I. INTRODUCTION 

The physics of systems driven through a random environment 
is by construction irreversible. The fluctuation dissipation re- 
lation does not hold and one expects the coarse grained de- 
scription to exhibit signatures of this irreversibility. In driven 
manifolds it has indeed been shown that non-linear Kardar- 
Parisi-Zhang (KPZ) terms are generated in the equation of 
motion, except when forbidden by symmetry [1,2]. A ques- 
tion which was debated for long time is whether at zero tem- 
perature these terms vanish as the velocity v — » + . This is 
the limit which is relevant to describe depinning (/ — * /+). 
It was found some time ago that there are two main universal- 
ity classes for interface depinning [3-5]. The conclusion was 
reached mainly on the basis of numerical simulations, which 
measure the interface velocity v(9) as a function of an aver- 
age imposed slope 9, as well as various arguments related to 
symmetry. In the first universality class, the isotropic depin- 
ning class (ID), the coefficient A of the KPZ term vanishes 
as v — ► + and the KPZ term is thus not needed in the field 
theoretic description. In the second class, the anisotropic de- 
pinning class (AD), v{6) still depends on 9 as / — > /+ and 
the KPZ term is present even at v — > + . For AD, numer- 
ical simulations based on cellular automaton models which 
are believed to be in the same universality class [6, 7], indi- 
cate a roughness exponent ( 0.63 in d = 1 and £ » 0.48 
in d = 2. On a phenomenological level it has been argued 
[6-8] that configurations at depinning can be mapped onto di- 
rected percolation in d = 1 + 1 dimensions, which yields 
indeed a roughness exponent Cdp = = 0.630 ±0.001, 

a dynamical exponent z = 1, a velocity exponent /?dp = 
mi — vj_ Rj 0.636 and a depinning correlation length exponent 



i^Dp = Mi = 1.733 ± 0.001. Some higher dimensional exten- 
sions of these arguments in terms of blocking surfaces have 
been proposed [9-12], but there is, to our knowledge, no sys- 
tematic field theoretical connection between these problems. 

Recently we have reexamined the functional renormaliza- 
tion group (FRG) approach, introduced previously [13-17] to 
describe isotropic depinning to one loop within a e = 4— d ex- 
pansion. We constructed [18, 19] a consistent renormalizable 
field theoretical description up to two loops, taking into ac- 
count the main important physical feature - and difficulty - 
of the problem, namely that the second cumulant A (it) of 
the random pinning force becomes non-analytic beyond the 
Larkin scale. The 2-loop result for the exponent £ shows de- 
viations from the conjecture[17] £ = (4 — d) /3. The reason is 
the appearance of "anomalous" corrections caused by the non- 
analytic renormalized disorder correlator. The 2-loop correc- 
tions proved to be crucial to reconcile theory and numerical 
simulations [18, 19]. 

The aim of this paper is to extend this FRG analysis to the 
universality class of anisotropic depinning. We first show that 
beyond the Larkin length, the KPZ-term is indeed generated 
at v = + , as long as it is not forbidden by symmetry. We 
explicitly compute the lowest order corrections for a simple 
model studied in recent simulations [20, 21]. Next we derive 
the FRG-flow equations for the second cumulant A(u) in a 
4— e expansion. In a previous study, Stepanow[22] considered 
the model to one loop, but did not take properly into account 
the non-analyticity of the renormalized disorder. Since this 
is physically important, we reexamine the problem here. In- 
deed, we find several new important "anomalous" corrections, 
including the one which generates the KPZ term in the first 
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place, as well as terms correcting the /3-function. We then in- 
troduce an equivalent description in terms of Cole-Hopf trans- 
formed fields. This description is not only much simpler to 
study in perturbation theory (e.g. to two loops it reduces the 
number of diagrams by an order of magnitude), but it allows 
us to obtain a number of results to all orders. We argue that 
the coefficient A/c which measures the strength of the KPZ 
non-linearity is uncorrected to all orders. We also determine a 
non-trivial subspace of the disorder correlators in the form of 
simple exponentials which is an exact invariant of the FRG to 
all orders. In the Cole-Hopf variables it is reformulated as the 
field theory of a specific branching process, or equivalently 
reaction-diffusion process. 

Our flow-equations allow to study both periodic disorder, 
relevant for charge density waves (CDW), and non-periodic 
disorder, relevant for lines or interfaces in a random environ- 
ment. In both cases we find several non-trivial fixed points. 
All these fixed points possess at least one unstable direction 
and should thus be associated to transitions. It seems that 
perturbatively the large scale behaviour is dominated by a 
runaway-flow, as it is in the standard KPZ problem[23, 24]. 
The difference is that its direction is a non-trivial function 
A(tt) in functional space. Analysis of the above mentioned 
invariant subspace suggests that the flow goes towards a spe- 
cific branching process. The present RG analysis is however 
unable to attain the non-perturbative fixed point. Thus, it also 
does not allow to strictly decide whether d — 4 is the up- 
per critical dimension of the anisotropic depinning problem, 
which is an open issue. 

Finally, since there are indications that KPZ terms may 
be needed in the description of the motion of a contact line 
[25], we have studied manifolds with long range elasticity and 
the simplest KPZ term. We determine the critical dimension 
above which this KPZ term is irrelevant, as well as the rough- 
ness at crossover. 

II. MODEL 

We consider a rf-dimensional interface (in d+ 1 embedding di- 
mensions) with no overhangs parameterized by a single com- 
ponent height field u(x). The case where the disorder is pe- 
riodic corresponds to a single component CDW in d dimen- 
sions. The common starting point is the equation of motion 

r]d t u xt = cd 2 x u xt + \{d x u xt f + F(x, u xt ) + f xt (II. 1) 

with friction 77, a driving force f xt — / and in the case of 
long-range elasticity we replace (in Fourier) q 2 u q by |<z| a Ug 
(with mostly a = 1) in the elastic force. The pinning force 
F(x,u) is chosen Gaussian with second cumulant 

F(x,u)F(x' >') = A(u - u)5 d {x - x) . (II.2) 

Temperature can be taken into account as an additional white 
noise rj(x,t) on the r.h.s. of (II. 1) with (r](x,t)ri(x' ,t')) = 
2rjTS(t - t')S(x - x'), but we will focus here on T = 0. 

Disorder averaged correlation functions (A[u xt \) = 
(A[u xt ))s an d response functions 5(A[u]) /6f xt = 



(u xt A[u]}s can be computed from the dynamical action 

S= u xt (r)d t - cdl)u xt - \u xt (d x u xt ) 2 (H.3) 

-~ / U xt U x t>A(u x t ~ U xt >) ~ / Uxt fxt 
1 Jxtt' Jxt 

The uniform driving force f xt — f > (beyond thresh- 
old at T — 0) may produce a velocity v = d t (u xt ) > 0, 
a situation which we study by going to the comoving frame 
(where (u x t) — 0) shifting u x t — > u x t + vt, resulting in 
f —* f — r) v - This is implied below: Each A is of the form 
&{uxt — Uxt' + v(t — t')), and we always consider the quasi- 
static limit v = + . Perturbation theory is performed both in 
KPZ and disorder terms, using the free response function 

(«,/«_ ? , t )o - R q , t -t> = rT'e-t*-''^ 2 /^ - t'). (0.4) 
III. GENERATION OF THE KPZ-TERM 

In this section we show how the irreversible (non-potential) 
KPZ term is generated, even in the limit v — > + , starting 
from a purely reversible equation of motion, where all forces 
are derivatives of a potential. 

Let us first consider the model recently studied numerically 
by Rosso and Krauth [20, 21], where the elastic energy is 
J x E(Vu x ), and e.g. E{6) = §8 2 + f-6» 4 . The relevant con- 
tinuum equation of motion is: 

f]d t u xt = E"(d x u x t)d 2 u x t + F(x,u xt + vt) + f-r]v (III.l) 

Note first that when C4 = 0, which corresponds to the 
isotropic depinning class with E{9) = |# 2 , the generation 
of the KPZ term is forbidden by the statistical tilt symmetry 
(STS), i.e. the invariance of the equation of motion under a 
shift Uxt — * Uxt + fx with f x = hx (or more generally the co- 
variance under an arbitrary f x ) [29]. When C4 ^ the model 
does not obey STS and the KPZ term is not forbidden, and 
indeed it is generated at finite velocity v > 0. This consider- 
ation alone is insufficient to show that it is still generated as 
v —> since in that limit the symmetry u —> —u should forbid 
it. Indeed, if one performs conventional perturbation theory 
with an analytic disorder correlator A(u), one does immedi- 
ately find that the KPZ term vanishes as jj — > + . However 
one needs a mechanism by which, as v — > + , the symmetry 
u — > —u remains broken. 

As we now show, this mechanism is provided by the non- 
analytic nature of the disorder. We know from studies of 
isotropic depinning [14, 16, 18, 19] that at T = the 
coarse grained disorder becomes non-analytic (NA) beyond 
the Larkin length [30]. We show below that this is also the 
case for the situation considered here. 

Using the techniques developed in Ref. [18, 19] the cor- 
responding perturbation theory, with a non-analytic A(w) be- 
comes (see figure III. 1 for notation) 
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FIG. III. 1 : The diagram generating the irreversible nonlinear KPZ 
term with one disorder vertex (notations are as in Ref. [18, 19]) and 
one C4 vertex (the bars denote spatial derivatives). 



= f /c-( t+t ') fc2 (fcV + 2(M 2 ) 

V Jt>oJt'>OJk 

xA'(u x , t+t , - u x . a + v(t + t')) , 

(IH.2) 

At T = 0, u x ,t has vanishing expectation value and the argu- 
ment of A' becomes v(t + 1'). Using that 

A(«) = A(0) + A'(0+)|«| + ^A"(0+)u 2 + ... (IH.3) 
A'(u) = sign(u) + A"(0+)u + . . . (IH.4) 

and observing that t, t' > 0, (III.2) can be written as 

5X = -^J JJ e ~ (t+t ' )k2 ( fcV + 2{kp ^ 

x (A'(0+) + A"(0+)v(t + 1') + 0{v 2 )) 

(IH.5) 

The leading term of this expansion, which is the only UV- 
diverging one for 4 > d > 2, is obtained by setting v = 0. 
Integrating over t, t' and using the radial symmetry in k gives 



5\ = -a 1 + 



d) L 



A'(0+) 
k 2 



+ 0(v) . (111.6) 



Similarly, there is a correction to c, which reads 



5c 



leading to 




5c = C4 ( 1 + 



k k 2 \ + d 



(III.7) 



(III.8) 



As will become clear below, the natural coupling for the KPZ- 
termisnot A, but the ratio A = A/c, which is corrected as [31]: 

SX = -C4 (l + (A'(0+) + AA(0)) ^ 1 (IH.9) 

Thus we have shown that the symmetry u — > —u which for- 
bids the KPZ term (e.g. in an analytic perturbation theory 
where A'(0) = 0), is broken here at v = + by the non- 
analytic term, and that a KPZ term is indeed generated at 
depinning. As in our previous study [18, 19] the only as- 
sumption is that the interface always advances forward (or 



that backward motion can be neglected in the steady state), 
supported in this single component model by no passing the- 
orems [16, 20, 21]. By providing a physical mechanism, this 
explicit calculation confirms the argument given in [4] based 
on a Larkin type estimate of the angle 9 dependence of the 
critical force. 

Note the sign of the generated KPZ term. Since A'(0 + ) 
is negative, A is positive as found in simulations [3, 4]. It is 
a bit counter-intuitive that the surface should become stiffer. 
Also it effectively corresponds to the generation of a positive 
average curvature. This is presumably through non-analytic 
coarse grained configurations of the string (in d = 1) since 
otherwise V 2 u = [Vu]q would grow as L which is un- 
physical, while cusps in u(x) allow for such a result. 

This model is only a particular case, which shows that the 
anisotropic depinning class is rather broad and not limited 
to anisotropic disorder. In general, unless they are excluded 
by symmetry, KPZ-terms will appear. One such case, corre- 
sponding to a flux line in 1 + 1 dimensions which moves per- 
pendicular to itself was considered in [4]. There disorder is 
anisotropic with correlators A x and A u for the pinning force. 
In the case of isotropic disorder A x = A u , exact rotational 
invariance (which in infinitesimal form reads u — > u + Ox, 
x — > x — 9u) should suffice to exclude the KPZ term. We 
have indeed checked this by adding to the above MSR-action 
with A = the non-linear terms of [4] 

SS = -J «*t [AV 2 u xt (Vu xt f + Bf(Vu 2 xt )] (IH.10) 

xt 

- j u xt u xt ' [C(Vu xt ) 2 +DVu xt Vu xt >] A(u xt -u xt i) ■ 



The generated KPZ term reads to lowest order 



5\ = 2(-A + C + 



(III. 11) 



Since the equation of motion of Ref. [4] for A x — A^ cor- 
responds to A = —1,D = C = 1/2, one checks to lowest 
order that the KPZ term is indeed not generated. Although we 
have not checked it further, it is clear that this property should 
extend to all orders. In the anisotropic class A, can a priori 
be of any sign. The argument given in [4] suggests that for 
the flux-line model A is positive when A^ < A h and negative 
for A x > Ah- Note that anisotropy by itself is not enough to 
generate the KPZ term, but that a non-linear and non-analytic 
disorder correlator is needed, and that this term will of course 
not be generated in a simple Larkin-type random force model, 
where A x and A^ are constants. 



IV. DIMENSIONAL FLORY ESTIMATES 

Before using analytical methods, let us indicate a simple 
Flory, or dimensional, argument which indicates how expo- 
nents for ID and AD can differ. In the absence of a KPZ term 
and setting u ~ the two static terms in the equation of 



motion scale as 



V 2 u~ x c ~ 2 



F(u x ,x) — x 



(IV. 1) 
(IV.2) 



Using F(u,x)F(u',x') <~ S(u - u')S d (x - x') for random 
field disorder gives the Imry-Ma value 



d 



(IV.3) 



which can be argued to be exact for the statics and is corrected 
by 0(e 2 ) terms at depinning. These types of arguments typi- 
cally give the exact result for LR correlated disorder, as the LR 
disorder part is not renormalized. It happens that this range is 
long enough for the statics but not for depinning; hence there 
is a correction at depinning which increases Note that it 
becomes again exact for depinning if the range of A in u or x 
is large enough (see e.g. the end of Section IV B in [26] and 
Appendix B). 

In presence of a KPZ term the latter scales as 



(V x u) 



2 2C- 
~ X s 



(IV.4) 



Supposing that it is relevant, it dominates over the elastic term. 
Balancing KPZ-term against disorder gives the modified Flory 
estimate 



4-d 



(IV.5) 



For d = 1 it yields (p = 0.6 versus £ = 0.63 observed in 
simulations [20], which is not bad an estimate for such a sim- 
ple argument. Again it is possible that if one increases the 
range of A the estimate (IV.5) becomes again exact, as is the 
case for standard KPZ (directed polymer) see Appendix B. 
Note however that it works with an upper critical dimension 
d = 4, which is an open question, and is thus merely indica- 
tive. 



V. FLOW-EQUATIONS IN PRESENCE OF A 
KPZ-TERM 

Let us start by deriving the FRG flow of A, c, i] and A to 
one loop starting from (II. 3). The KPZ and disorder terms are 
both marginal in d = 4 and become relevant below. Simple 
dimensional arguments show that these are the only needed 
counter- terms. We have computed the effective action to low- 
est order. The corrections as given by the diagrams on figure 
VI are (for details see Appendix A): 



Si] 

V 
5c 

c 

5X 



[a c- 3 AA'(0+) + c- 2 A"(0+)] / 
= - [aiAc- 3 A'(0+) + a 2 A 2 c- 4 A(0)] / 



(VI) 



— = - [a 3 Ac- 3 A'(0+) + a 4 A 2 c- 4 A(0)] / 
A 

SA = [a 5 A 2 c- 4 A 2 + c- 2 (A"(A(0) - A) - (A') 2 )] I 




FIG. V.l: 1 -loop-diagrams correcting A (top), c (middle), and r\ and 
A(w) (bottom). 

where I = J 1/q 4 (integrated over the shell if using Wilson's 
scheme) and the coefficients are: 



a = 1 , ai 
as = a,4 = 4/d , 



2(d-2)/d, 
a 5 = 2 . 



a 2 = 4/d 



(V.2) 



In the following we will set d = 4 in these coefficients since 
they are universal only to this order. This gives 



d0 = a l = a 2 = «3 = &4 = 1 



a 5 = 2. 



(V.3) 



One then notes that the quantity A/c remains uncorrected to 
first order in d = 4. In the next section we shall argue that this 
remains true to all orders. The corrections to the linear term in 
(II. 3) can be interpreted as the correction to the critical force: 

Sf = -6f c = (Ac- 2 A(0) + c- 1 A'(0+))/ 1 (V.4) 

where I\ = J q ^ . It does not require an additional counter- 
term if we tune / to be exactly at depinning / = f c . 

In view of the non-renormalization of A/c in (V. 1) it is use- 
ful to denote the unrescaled coupling constants as 



A = - 

c 



T}=- 



(V.5) 



One should also notice that if one performs the change of vari- 
able in the initial model u — > u/X, u — > uX, then the free 
(quadratic) part of the action (proportional to c and rj) remains 
invariant while disorder and KPZ terms become: 



A-»l 
A(u) -» A 2 A(w/A) 



(V.6) 



Thus the coefficient A can be set to one upon appropriate re- 
definitions of disorder and displacements. 

It is natural to start the study of the FRG flow and the search 
for fixed points as for A = by defining the following rescaled 
parameters 



A 

A(u) 



AA7< 



:Af- e A(uA7<) 



(V.7) 
(V.8) 
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within a Wilson scheme where = Ac~ £ is the running UV 
cutoff. This yields two coupled equations for the couplings A 
and A(u) 

<%lnA = C (V.9) 
d e A(u) = (e - 2C)A(u) + «CA'(u) 

+2A 2 A(u) 2 + [2A 2 A(0) + 2AA'(0+)]A(u) 
-A'(u) 2 - A"(u)(A(u)- A(0)) (V.10) 

where here and below we absorb el — S4/ (2ir) 4 in the cou- 
plings. One notes that if there is a fixed point for A(u), then 
£ is the roughness exponent since 



2„4 



(u q u_ q ) = A{Q)/c 2 q 



= Ap*A*(0)/<T~ A*(0)/<?' 



<i-2C 



(V.ll) 



when evaluated at scale Kg = q. A more rigorous calculation 
uses the effective action[19] at non-zero momentum, but to 
one loop gives the same result. The dynamical exponent z in 
t ~ x z and the anomalous dimension of the elasticity can be 
determined from 



-iP = d e \nc= -AA'(0+) - A 2 A(0) 
- 2 = d e lnfa/c) = -A"(0+) + A 2 A(0) 



(V.12) 



The correlation-length exponent v in £ <~ (/ — f c )~ v and the 
velocity exponent (3 in v ~ (/ — f c )^ are given by the scaling 
relations 



2-C + ^ 



2-C + V 



(V.13) 
(V.14) 



This can be seen by noting that the action (II. 3) is invariant 
under x — e e x', t — e ze t', u — e^u', u = { l ' c ( 2 - z -C- d +'>P) e 
provided 77 = rj' e ( 2 - z +<W, c = c ' e <^, A = A' e^ +c ^, 
/ = f'c^-t+W and A = A' e ( £ - 2 C+ 2 ^ as well as T = 
T / c (2-d-2C+</>)£ while in presence of STS one has ip = 0, 
this is not the case here. In a Wilson formulation, the critical 
force is obtained by integration over scales of 



5£/ c = -q(A £ A,(0) + A^(0+))A 2 



2-C 



(V.15) 



a quantity which physically is likely to remain positive. 

A salient feature of the AP class is that the critical force 
depends on the angle by which the interface is tilted. From the 
arguments of [3, 4] the characteristic slope 9 should scale like 
the ratio of the characteristic lengths orthogonal and parallel to 
the interface, 9 ~ ~ (/— fc) v ^~^ and more generally 

the velocity should behave as 



v(f,0) = (f - fc(0)f 9 



. (V.16) 



Defining A e ff by [4] v(f, 9) = A e ff6> 2 + . . ., the small 9 expan- 
sion of v(f, 9) gives he effective A ff as 

Acs ~ (/ - fMf~ 2v{1 ~ i] = (/ - W)r v{2 - C - Z) ■ 

(V.17) 



Performing the redefinition u = u + 6x, we can compute the 
critical force as a function of the angle to lowest order in 
disorder 



6fc{8) 



= -9 2 X 1 + 



d 

T 



(a 2 A(0) + AA'(0+)) 



(V.18) 



and thus we find an angular dependence, which is increased 
under renormalization. 

The notable feature of the above FRG equation is the ab- 
sence of corrections to A to this order in eqs. (V.l). It is cru- 
cial to determine whether this persists beyond one loop. If 
there were corrections to higher order this might allow for a 
non-trivial fixed point of A and thus to fix £. On the other 
hand, absence of corrections would imply that for ( > 0, A 
flows to infinity, which makes the existence of a perturbative 
fixed point doubtful. In the next section, we present a different 
approach, which allows to clarify this question. 

It is worth noting, that since KPZ-terms are only gener- 
ated above the Larkin length, the FRG flow below the Larkin 
length (as well as the value of this length) is identical to the 
case A = 0. It is however instructive to artificially consider 
the above FRG flow for an analytic function and with a given 
imposed bare value of A (setting ( — 0). One gets 

di A(0) = eA(0) + 4A 2 A(0) 2 (V. 19) 

d e A"(0) = eA"(0) - 3A"(0) 2 + 6A 2 A(0)A"(0) . (V.20) 

The bare disorder has A(0) > and A"(0) < 0. Since all 
terms on the r.h.s. of (V.20) have the same sign, |A"(0)| di- 
verges faster if A 7^ 0, meaning that the KPZ-term cannot 
prevent A(u) from becoming non-analytic. Note that the first 
equation exhibits a runaway at -£<a(o) which can shorten the 
Larkin length. Ind = 4 + eatA = there is an unstable fixed 
point at A"(0) = — e/3 separating a Gaussian weak-disorder 
phase with the bare unrescaled Larkin force producing finite 
displacements, and a phase where disorder seems to become 
non-analytic, only to become irrelevant at larger scales as can 
be seen by examining the flow in the non-analytic space be- 
yond the Larkin length. At A > there is a fixed line at 
A(0) = — e/ (4A 2 ) > which separates a phase where A(0) 
grows from a phase where it decays to zero. On the transition 
line the flow is towards a non-analytic disorder. 

VI. COLE-HOPF TRANSFORMED 
THEORY 

We now introduce the Cole-Hopf transformed theory which 
has a lot of interesting properties. 

Starting from (II. 1) we first divide by c. This gives 

fjd t u xt = d x u xt + X{d x u xt ) 2 + -F(x, u xt ) + - (VI. 1) 

c c 

We then define the Cole-Hopf transformed fields 

Z xt :=e^ * Uxt = X ^A. (VI.2) 

A 
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The equation of motion becomes after multiplying with \Z xt 
rjd t Z xt = d 2 x Z xt + —F (x, Mf^l ) Zxt + ^-Z xt (VI.3) 



and the dynamical action 

S = [ Z xt (ijdt - d 2 x ) Z 

Jxt 

I Z xt Z xt Al 

1 Jxtt' \ 

c J x 



xt 

In Z xt - In Z xV 
A 



Z X ±Z. 



xt^xt 



Zxt'Z x t' 
(VI.4) 



It is important to note that the above formal manipulations 
are only valid in the mid-point (Stratonovich) discretization. 
The strategy therefore is to start from the original equation of 
motion, which is interpreted in the Ito discretization, switch to 
Stratonovich, make the change of variables, and then switch 
back to Ito. Note the identification: 



Uxt = —Z x tZ x t 
c 



(VI.5) 



and that in this formalism the force (or the distance to the 
critical force) corresponds to a mass: 



m 2 = -(/-/ c ) 
c 



(VI.6) 



Let us first illustrate how perturbation theory works in this 
new formulation and how one can easily recover the 1-loop 
FRG equation obtained in the previous section. Perturba- 
tion theory is performed with the standard response-function. 
We note a very important property: To contract Zqq with a 
disorder-insertion Z xt Z xt \ 2 A^ lnZxt ^ lnZxt ' ^ Z xt iZ xt i and 
focusing on Z xt (not Z xt >), one can decide to either contract 
Z xt standing outside the A or inside. In the first place, this 
eliminates the factor Z xt , but leaves A underived. In the sec- 
ond case, deriving the argument of A, gives A'/ A, together 
with a factor of 1 jZ xt from the inner derivative. The latter 
also cancels the Z xt standing outside the A. So independently 
of where one derives, one always looses the factor of Z xt out- 
side A. Contracting n times towards the vertex at x, t thus 
gives a factor of Z x ^ n . This observation shows that the dia- 
grammatics are a very simple generalization of the case with- 
out the KPZ-term which was detailed up to two loops in [19]. 
One easily verifies that the latter case is reproduced upon con- 
tracting only the argument of A. To see this, one performs the 
perturbation theory and finally takes the limit of A — > 0. Each 
time, one has contracted a Z xt outside of A, one is missing 
a factor of 1/ A, and the term vanishes in the limit of A — > 0. 
Further remark that for A — > 0, the argument of A becomes 



Z x t — Z x f 
A 



u xt - u xt > + 0(A) 



(VI.7) 



This shows that the perturbation theory for isotropic depin- 
ning is reproduced. 



Thus the new diagrams, in the presence of the KPZ-term, 
can be deduced from those for A = by allowing additional 
contractions of a Z xt outside the A. Compared to performing 
calculations using (II. 3) this yields a much simpler perturba- 
tion theory, with far less distinct diagrams. E.g. to two loops, 
the number of diagrams is reduced by at least a factor of ten. 

Note that now a renormalization of the term ZKZ is al- 
lowed, since it is no longer forbidden by STS. Indeed shifting 
Uxt — > u xt + ax/X and Z xt — > Z xt e~ ax , we find that the 
action changes by 



5S 



= [ Z xt (a 2 + aV) Z xt . 

Jxt 



(VI.8) 



However, since the action (VI.4) is still translationally invari- 
ant, it remains unchanged under 



Z x t 
Z x t 



LiZ xt 

x z 



(VI.9) 



Transforming only Z xt — » nZ xt without changing Z xt will 
allows us later to fix the coefficient of the Laplacian to unity 
and transfer all its corrections into corrections to A and fj. 

We now present the calculations at 1 -loop order. We start 
with the corrections to fj. Contracting one disorder vertex 
once with itself, we obtain 



A Z x t'Z x 



In Z xf i— In Z r 



A 



+ iA'f ,uZ " 



A 

Expanding In Z xt > — In Z xt for small times yields 



lnZ r 



(VI. 10) 



In Z xt , - In Z xt = {t ' f J dtZxt + Q(t - t'f (VI.11) 



Jxt 



One also has to expand Z xt around xt': 

Z x t = — (t — t)d t Z xt ' 



(VI. 12) 



Since the manifold only jumps ahead, the arguments of A and 
A' are always positive. Putting all terms together, we obtain: 



Z x t'dtZ x t'(t' — t)Ro tl 



/ AA'(0+) + A"(0+)) - (A 2 A(0) + AA'(0+)) 
Integrating over t' — t yields 



Z x t'dtZ x t'I 



(AA'(0 + ) + A"(0+)) - (A 2 A(0) + AA'(0+)) 



(VI.13) 



(VI. 14) 



We have grouped terms such that in the first bracket there ap- 
pear the corrections to — ^ and in the second those to ^. 
Here they appear all together in one diagram. In the absence of 
the KPZ-term only the term independent of A survives. Not- 
ing the cancellation between the two terms, we finally arrive 
at 



Sfj 
V 



A"(0+) - A 2 A(0) 



(VI. 15) 
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We now turn to corrections to disorder. Reminding that the 
arrows can either enter into the argument of A or into the sin- 
gle Z-field, we get the following contributions (plus some odd 
terms, which we do not write): 



5A(u) a = 
SA(u) h = 
SA(u) c = 
5A{u) d = 



-A"{u)A(u) + X 2 A(uf 



-A» 2 + X 2 A(uf 



A"(u)A(0) 



(VI. 16) 



I 



AA(w)A'(0+) + A 2 A(u)A(0) 



I 



These reproduce the corrections obtained in the previous sec- 
tion, but quite differently. 

The Cole-Hopf transformed theory suggests that 



sx = o 



(VI. 17) 



to all orders. To prove this one has to show that the following 
terms are not generated in the effective action 



z x t-^—(yz x t) 2 . 



(VI. 18) 



It is easy to see that these terms result from a change of A 
(keeping u xt and u xt fixed): 



Z x t — * Z xt ^1 + -7T- In Z xt 
Z x t — * Z xt ( 1 j- In Z xt 



(VI. 19) 
(VI.20) 



and thus the Laplacian generates (VI. 18). One can also again 
consider a term like C4 which is known to produce a shift in 
A (see (III.9)), and does produce (VI. 18) above together with 
other irrelevant terms with more gradients. In fact (VI. 18) is 
by power counting the only term marginal in d = 4 which can 
appear. This term could in principle come from vertices with 
several derivatives acting on A at point x. As previously dis- 
cussed, it is always compensated, but the compensating factor 
could be on a different vertex at position x' and hence produce 
(VI. 18) via a gradient expansion. We have shown in Fig. VI. 2 
the 2-loop diagrams correcting terms with a single response 
field in the effective action and the Z and 1 jZ fields which 
appear at each vertex. All terms contribute to fj. Graphs b, c 
and d each give a term of the form (VI. 18) by expanding the 
Z 2 on the lower disorder, but the sum of them cancels. As 
we will discuss below this is graphically achieved by moving 
the ends of the arrows around on the upper vertex, suggesting 






z* 




f z 





FIG. VI. 2: 2-loop dynamical diagrams correcting the single Z- 
component. Diagrams a - g correct the friction. Only diagrams 
e and f have a sufficiently strong divergence in space (after time- 
integration) that they can produce spatial gradients. In fact they both 
correct Z&.Z. (The diagram is the well-known sun-set diagram from 
<?i 4 -theory.) 






FIG. VI. 1: 1-loop dynamical diagrams correcting A 



FIG. VI. 3: Figure explaining the non-renormalization of A, see main 
text. 



a more general cancelation. Another argument is that the di- 
vergence in space between the upper and lower vertex is not 
strong enough in order to contribute to (VI. 1 8) or J ZAZ. For 
this to happen, one needs three response-functions between 
upper and lower disorder, as is the case for diagrams e and 
f. They thus both contribute to J ZAZ, but since they have 
only a single Z on the lower disorder, they do not contribute 
to (VI. 18). 

We now argue that to all orders in perturbation theory no di- 
agram proportional to a single Z (one connected component) 
can be generated, which contains a factor of (VZ) 2 ^. We 
believe these arguments to be conclusive; especially we have 
not been able to construct any counter-example at 3- or4-loop 
order. However the structure of the theory is sufficiently com- 
plicated that some caution is advised. 

Look at figure VI. 3. The response-functions (arrows) in an 
arbitrary diagram correcting a single-time vertex have a tree- 
structure (left). This diagram can be completed by adding the 
disorder-interactions between arbitrary pairs of points (mid- 
dle). A potentially dangerous factor of appears at point 2. 
Point 2 has a "brother" 3, to which it is connected by a disor- 
der correlator A (dashed line). 

Then, two cases have to be distinguished: Either there is 
no line entering point 3, then point 3 can contribute his factor 
of Z to point 2: Since it is at the same point in space, the 
difference can be expanded in a series in time, giving time- 
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derivatives of Z which do not spoil the argument. 

On the other hand, there may be a line entering point 3. 
This is drawn on figure VI. 3 (middle). By construction (at 
least) two branches (of response-functions) enter at point 2. 
At least one of them does not contain the brother of 2 (here 
point 3). Here it is the left branch, containing point 1. Now 
consider the diagram where the response-function from 1 to 2 
is replaced by a response-function from 1 to 3 (right). Since 
one can always contract last the response-field at point 1, lead- 
ing to either the response-function from 1 to 2 or the one 
from 1 to 3, these diagrams have the same combinatorial fac- 
tor, but differ by a a factor of —1, due to the derivative of 
A( lnZ ~^ nZ ) on either the first or the second argument. This 

comes in both cases with a factor of ^, at the same position in 
space but at different positions in time. However, due to the 
tree-structure, the time-integration can always be done freely, 
and the two vertices finally cancel. This argument is sufficient 
before reaching the Larkin-length. However after reaching 
the Larking length, the non-analyticity of the disorder may 
yield additional sign-functions in time between both ends of 
the vertex, as has been observed in [19]. Then the proof gets 
more involved. There is another very powerfull constraint on 
the generation of terms like (VI. 18): One has to construct 
a diagram with a strong spatial ultraviolet divergence, such 
that after Taylor-expanding Z in space the additional factor 
of x 2 together with this strong ultraviolet divergence gives a 
pole in 1/e, i.e. a logarithmic divergence at d = 4. This is 
the situation for diagrams e and f in figure VI. 2. It arises if 
and only if there are 2n + 1 response-functions connecting n 
points in space (this may well be a sub-diagram), but where 
response-functions that connect the same point in space are 
not counted. In all examples which we considered up to 4- 
loop order, which had sufficiently many factors of 1 /Z, and 
which had the correct UV-structure, the (2n + 1) response- 
functions where enough to enforce an ordering of times, such 
that the mounting proof sketched on figure VI. 3 went through. 
We have to leave it as a challenge to the reader to either find a 
counter-example or to make the above arguments rigorous. 

Let us now return to the analysis of the RG-equations. We 
introduce rescaled variables according to 

A(u) = A^ 2C A(uA<;) (VI.21) 
A = AA7 C (VI.22) 

with A? = Ac~ e . Because we have defined Z = c Xu , in order 
not to generate additional terms, a rescaling of u demands a 
(compensating) rescaling of A such that the product remains 
unchanged. Even though this may not be the best choice cor- 
responding to the existence of a fixed point, it is the only way 
to preserve the Cole-Hopf-transformation, leaving Z and In Z 
unchanged. The rescaling of A comes from the rescaling of 
A, which appears as a factor of A 2 in front of A in the action 
and as a factor of 1/ A in the argument of A. 

This leads again to the FRG flow equation given in (V.10): 

d t A{u) = (e - 2C)A(u) + C«A'(«) 

-A"(u) (A(u) - A(0)) - A» 2 



+2AA(u)A'(0+) 

+2A 2 (h{uf + A(u)A(0)) (VI.23) 

Further remarkable properties of the Cole-Hopf transformed 
theory will be shown below. We now turn to the study of the 
FRG flow. 

VII. PERIODIC CASE 

We now consider the case, where A(u) is a periodic function 
with period 1. The starting point is (VI.23) with £ = 0, thus 
A = A remains constant under renormalization (to all orders). 
Since the period is fixed, A cannot be scaled away using (V.6). 
It is thus a continuously varying parameter and we must study 
the flow as a function of it. 

In eq. (VI.23) there is a tendency for a runaway flow, as 
can be seen by analyzing the flow-equation (VI.23) with the 
trivial solution A(u) = A 

d e A = eA + 4A 2 A 2 . (VII. 1) 

This corresponds to the localization - or self attracting chain 
- problem studied in [27] and we expect on physical grounds 
the full functional form of A(u) to be important, which may 
lead to other fixed points. 

For A = we already know that there is an unstable fixed 
point 

A e (u) = A* (u) + c e ee (VII.2) 
A*(u) = ^-^u(l-u), (VH.3) 

which describes isotropic depinning for CDW. This fixed 
point survives for small A as can be seen from a series ex- 
pansion in powers of A. Moreover at each order in A, A*(u) 
remains polynomial in u(l — u). We do not reproduce this ex- 
pansion here, since we have succeeded in obtaining the fixed 
point analytically. Equation (VI.23) possesses the following 
remarkable property: 

A three parameter subspace of exponential functions 
forms an exactly invariant subspace. 

Even more strikingly, this is true to all orders in perturbation 
theory. This property, which is quite non-trivial, is understood 
in the Cole Hopf theory, as discussed below. 
For our purposes, it is more convenient to write 

A(u) = ^-ef(u\) , (VII.4) 
A^ 

such that / satisfies the same FRG equation (VI.23) with A = 
e = 1, but with period A. This allows to make an ansatz for a 
family of exponential functions 

f(u) = a + be- u + ce u . (VH.5) 

The FRG-flow (VI.23) closes in this subspace, leading to the 



A = 2 



A = co 



FIG. VII. 1: Fixed point structure for different values of A. The coordinate system is such that a grows to the right and b to the top. The both 
separatrices are b — — ^ae A (blue/dark) and b — —a/(l + e~ A ) (orange/bright). 



simpler 3 -dimensional flow: 



dia = a + 4a 2 + 4ac 
d e b = b{l + 6a + b + 
dec = c(l + 6a + b + 



>rAbc 

5c) 

5c) 



(vn.6) 
(vn.7) 

(VII.8) 



This works only for amplitude one in the exponential; other- 
wise higher modes are generated. Also note that these equa- 
tions are not symmetric under the exchange of b and c, as one 
might expect from the interpretation we will present later. 

Requiring periodicity, or equivalently f(u) — f(X — u) 
imposes 



be 



(VII.9) 



and one checks that b/c is indeed unrenormalized. Thus one 
can study the simpler 2-dimensional flow 

x (VII. 10) 
(VII. 11) 



dg a ■ 

dib- 



a + 4a 2 + Aabe 
Ml + 6a + b + 5be 



A + 4&V 



as a function of A. A physical requirement is that 

A(0) =a + b(l + e~ A ) > (VII. 12) 



/ c ~-(A'(0) + AA(0)) 

= -\(a + 2bc- x ) > 
A 



For a < this is possible only if 

a 



1 



-A 



<b< --c x 



(VII. 13) 



(VII. 14) 



On the other hand, for a > the flow for a is always a — > oo 
in a finite time. Indeed the r.h.s. of (VII. 10) is always positive 
for a > 0. For b > this is trivial; for b < this can be seen 
from 



a + 4a 2 + 4a6c 



— A 



46 2 e" A 

^2 „-A 



= a + 4 (a + by c~ A - 4a6e~ A + 4a 2 (l - c~ A ) 
>0 (VII. 15) 

The flow given in (VII. 10) and (VII. 11) is shown in figure 
VII. 1 . There are four fixed points for d = 4 — e. In the original 
variables they are 

(i) Gaussian fixed point G (repulsive in all directions) with 

A(u) = 0. 
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(ii) Self-avoiding polymer fixed point SAP, where the corre- 
lator is a negative constant: 



A(u) 



e 
4A 2 



(VII. 16) 



It is the problem of localization in an imaginary random po- 
tential, i.e. the Edwards version of the better known self- 
avoiding polymer. It is attractive in all directions, even those 
not drawn here. Writing f(u) = —l/A + <p(u) and linearizing 
(VI.23) gives 



cw(o) = -m 



^'(0+) 



(VII. 17) 



d t <l/(u) = --<(>' (u) 



(VII. 18) 



This self-avoiding polymer fixed point will not play a role in 
the following since for the disordered problem A(0) > 0. 
However it is interesting in other contexts, as discussed below. 



(iii) Fixed point U, with one attractive and one repulsive direction. 



A(u) 



A 2 



1 + 54e- A + 5e- 2A - (1 + 5e- A ) yjl + C - A (34 + q- a ) 



+ 



8(l-5c- A (c- A -8))) 
2 



(1 - 7c- A - 3-v/l + c- A (34 + e" A ) 



( e -Au + e -A(l-ti)) 



The value at zero 



A(0) 



3 + e A (3 + Vl + c- 2A + 34c- A ) 
2A 2 (7 + c A (3Vl + c- 2A + 34^ - 1)) 



(VII. 19) 



(VII.20) 



is always negative for A > 0, thus the FP is unphysical for our problem in d = 4 — e. The combination yielding the corrections 
to the critical force 



fc = 



-1 + c A (7 + Vl + c~ 2A + 34c- A ) 
2A(7 + c A (3Vl + c- 2A + 34c- A - 1)) 



(VII.21) 



is always positive for A > 0. 

(iv) The random periodic fixed point RP has: 



A(«) = 1 



1 + 54e~ A + 5c~ 2A + (1 + 5c- A ) v /l + c- A (34 + c- A ) 



8(l-5e^ A (c- A -8))) 



(1 - 7c- A + 3-v/l + c- A (34 + e" A ) 



(c- Au + e 



(VII.22) 



A(0) = 



3-e A (-3 + VT 



34e 



2A 2 (-7 + e A (l + 3Vl + c- 2A + 34c- A )) 

(VH.23) 



f c ~ -(A'(0+) + AA(0)) 

_ -7 + c A (l + Vl + c- 2A + 34c- A ) 
_ 2A(-7 + c A (l + 3Vl + c- 2A + 34c- A )) 



(VII.24) 



Both quantities A(0) and f c are positive for all A > 0, thus 
this fixed point is physical. 

The fixed point RP is the continuation of the fixed point 
(VII.3) at A = 0: Note that apart from a constant only the 
term u(l — u) survives from the exponential functions. Like 
the fixed point at A = 0, it is attractive in one direction (to- 



wards the fixed point SAP) and repulsive in another (towards 
large A(m)). It is thus a critical fixed point. One can argue 
that any perturbation which leads to SAP is unphysical, since 
at some scale A(0) becomes negative. Since we did not find 
any strong reason why the system would be exactly on this 
critical surface, it is more likely that this FP represents a criti- 
cal regime which lies on the boundary of the physical domain. 
It is however interesting that its analytic form can be obtained. 
In particular one can compute correlation functions exactly at 
RP. 

An important question is whether there are fixed points out- 
side of the exponential subspace considered above. Let us give 
a few general properties. First the flow equations and fixed 
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point conditions near u = 

d e A(0) = C A(0) + 4A 2 A(0) 2 - A'(0) 2 + 2AA(0)A'(0+) 
*A'(0 + ) = A'(0+)(e + 2AA'(0+) + 6A 2 A(0) - 3A"(0+)) 

(VII.25) 

and the flow equation for J A 

d e J duA(u) = ( 



e + 2AA'(0+) + 2A 2 A(0) 
1 

+2A 2 / duA(?i) 2 



du A(u) 



(VII.26) 



shows that starting from J A = 0, a positive value for / A is 
generated in the early stage of the RG. If there is a fixed point 
value for J A it must be equal to 



Jo 



duA*(u) 



2A 2 J duA*(u) 2 



e + 2AA*'(0+) + 2A 2 A*(0) 



(VII.27) 



For small A at least this appears to be negative and of 0(A 2 ). 
From the flow-equation for A' (u) 

d e A'{u) = -£'"(«) (A(u) - A(0)) 
+A'(u) \e + 2AA'(0+) + 2A 2 A(0) + 4A 2 A(u) - 3A"(u) 

one sees that the behaviour at large u/X must be exponential. 
It seems that there are no non-exponential fixed points. 
The runaway flow will be discussed in the next section. 

VIII. RANDOM FIELD DISORDER 

Let us now consider non-periodic functions. The main prob- 
lem with the natural rescaling of u = u'e^ as in (VI. 22) is that 
A grows exponentially, and no fixed point can be found. Let 
us therefore study (VH.6)-(VH.8) setting the rescaling factor 
C = 0. Again we consider the invariant subspace of exponen- 
tial functions, parameterized by 



A(u) = ^e/(uA) 
f = a + be~ u 



(VIII. 1) 

(vin.2) 



for u > 0. Note that we have put the coefficient c = 0, since 
we are not interested in solutions growing exponentially in u. 
The flow is 



da = a + 4a 2 

db = b(l + 6a + b) . 

The physical requirements now read 

A(0) ~a + b > 
f c ~ -a > . 

So it is natural to look in the regime 

b > -a > 



(VIII.3) 
(VIII.4) 



(VIII.5) 
(VIII.6) 



(VIII.7) 



There is again the fixed point 

/(*) = - 



1 1 



(vni.8) 



which is the infinite A limit of the fixed point RP of the pre- 
vious section. Since f(x) does not go to zero at infinity as 
is expected for random field disorder, and since it is unstable 
along the line a = —1/4 it is unlikely to have any physical 
relevance for the anisotropic depinning class. The other fixed 
point is 



/(*) = -e" 



(VIII.9) 



which has the wrong sign. One clearly has runaway-flows 
within the exponential subspace. 

We have examined the flow of the FRG numerically. For 
all initial conditions considered, which were not exactly at 
one of the fixed points mentioned above, we found the so- 
lution to explode at some finite scale, a phenomenon which 
is known as the Landau pole. One issue is to identify the 
corresponding direction in functional space. This issue is re- 
lated to fixed points in d = 4 + e dimensions which we now 
briefly address. The diagram for 4 + e is obtained by changing 
A — > —A and de — ► — d^. This means to replace a — > —a and 
b — > —6 on figure VII. 1 as well as inverting the direction of 
all arrows. U then controls the boundary between the strong 
coupling regime of KPZ and the Gaussian fixed point G; SAP 
between localization (attractive polymers); the Gaussian fixed 
point is multi-critical and RP between branched polymers and 
Gaussian. For the random field case we now have 

A(«) = ^-ef(uX) (VIII. 10) 



The fixed point RP gives 



and the fixed point U is 



/(*) - e" 



(VIII. 11) 



(VIII. 12) 



which has the correct sign. It has a vanishing critical force, 
but is a good candidate for the critical behaviour between the 
Gaussian phase and the strong coupling KPZ phase. 

Let us now study the runaway flow for d = 4 — e. Suppose 
that Ai(w) is the solution of the (4 + e) -dimensional flow 
equation at e = 1. Then 

A/(u):=<fcAi(u) (VIH.13) 

leads to the flow-equation for the amplitude gi 

d gi = € gi + g 2 e . (VIII. 14) 

For the RF-case one has one such point at the boundary of the 
physical domain, as can be seen from the flow-equations 



a = Q 
db = b + b 2 



(VIII. 15) 
(VIII. 16) 
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Also not that since this mode explodes after a finite renormal- 
ization time, it is difficult to avoid. However, we have not yet 
completely ruled out another scenario, where at least some 
trajectories have exponential growth. Making the ansatz 



A l (u)=e(e 2 «f(u)+g(u)), 



(VIII. 17) 



this requires to find a solution to the /3-function at e = 0, 
which we write symbolically 



/?(/,/) = o. 



(VIII. 18) 



One can check that near zero such a solution is in principle 
possible. There is a solution, which vanishes at u = u* = 
1.39895 (for A = 1) and becomes negative beyond. One can 
argue that one needs it only up to u = uo < u*, since the 
linear term can no longer be neglected when f(u) approaches 
zero. Notingr = (l + V5)/4onehas /'(0) = -l/r,/"(0) = 
( 1 + 2r 2 ) / (3r 2 ) . In this scenario ( is determined together with 
g. It is unclear how this carries to higher orders, since it seems 
to require that f(u) is also solution of the /3-function at e = 0. 
This is however exactly what happens in the case A = with 
the constant shift A(0). Although numerics does not seem to 
confirm it, it is hard to disprove. A question which remains to 
be answered is what the basin of attraction of runaway growth 
and eventually of exponential growth are. 

IX. GENERAL ARGUMENTS FROM THE 
COLE-HOPF REPRESENTATION AND 
BRANCHING PROCESSES 

In the Cole-Hopf representation, it is easy to see why the ex- 
ponential manifold is preserved to all orders. Let us insert 



A(u) = -^(a + bc- Xu + ce Xu ) 



(IX. 1) 



in (VI. 4). The complicated functional disorder takes a very 
simple polynomial form 



S= Z xt (ijdt - d 2 x ) Z. 



X I 



X Jt<t' 



(aZ xt Z xt ,+bZ 2 t + cZ 2 t ,) .(IX.2) 



Note that we have ordered the vertices in time to distinguish 
between b and c taking correctly into account that the full cor- 
relator for the present non-analytic e.g. random field, problem 



time 



FIG. IX. 1: The three vertices proportional to a, b and c in 

fx t<t> Z*tZ xt i (aZ xt Z xt , + bZ xt + cZ xt ,). 



FIG. IX.2: Diagrams correcting the disorder in the branching 
representation. 



is (IX. 1) with u replaced by |u| (if (IX. 1 holded as an analytic 
function there would be no distinction between b and c, "thus 
no arrow of time"). 

The vertices presented on figure IX. 1 can be interpreted as 
branching processes, and we shall thus call this form branch- 
ing representation. Let us show how one reproduces the flow- 
equations (VII.6)- (VII. 8). In the time-ordered representation, 
diagrams a to d of figure VI. 1 have the form given on figure 
IX.2. To simplify notations, we set A = 1. Then 



A(u) =a + bc~ u + ce u 
A' (u) = -bc- u + ce u 
A" (u) = be~ u + ce u 

A(0) = a + b + c 
A'(0+)=c-6 
A"(0+)=6 + c 

The diagrams have the following contributions 



5A a (u) 



(1X3) 
(IX.4) 
(IX.5) 
(IX.6) 
(IX.7) 
(IX.8) 



5a = 


aal 






Sb = 


abl 




(IX.9) 


Sc = 


acl 






5a = 


(aa + 46c)/ 






Sb = 


2bal 




(IX. 10) 


5c = 


2acl 






5a = 









5b = 


26(a + 6 + c)| 


(from c 2 ) 


(IX. 11) 


5c = 


2c(a + 6 + c)| 


(from Ci) 




5a = 


2a(a + 2b) 


(from di and d 2 ) 


5b = 


2b(a + b + c)I 


(from d 2 ) 


(IX. 12) 


5c = 


2c(a + b + c)I 


(from di) 





5A h (u) 



5A c {u) 



5A c (u) 



Note that the factors of 2 come in general from contracting 
Z xt . The non-trivial factor of \ is due to the fact that the two 
right-most points in ci and c 2 are time-ordered. To relate the 
integral to /, one can first symmetrize (yielding the factor of 
i) and then freely integrate over time. Also note that only the 
last diagram, di + d 2 contributes to the asymmetry between b 
and c. 

In the same way, one can reproduce the corrections to 77. 
The only vertex in (IX.2) which contributes at leading order 
is the one proportional to a: b does not allow for a contraction 
and c will have both Z and Z at the same point, thus only 
corrects the critical force, a leads to 



Z x tZ x t'R x ,t-t> 



(IX. 13) 



and after a gradient-expansion following the procedure de- 
scribed after (VI. 10) we have 



Z x tZ x t 



(IX. 14) 
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FIG. IX. 3: The three phases of the flow-diagrams on figure VII. 1. 

Integration over t' leads to the correction to fj 

^- = -aI, (IX. 15) 

V 

which is the same one obtained from (VI. 15) using (IX. 6) and 
(IX.8). 

Let us now exploit this representation further: It is immedi- 
ately clear, that one cannot generate e _2A " which corresponds 
to 

/ / Z xt Z xt & (IX. 16) 

Jx Jt<t' ^xt 

or any other such fractions. This shows that the space of 
functions spanned by (IX. 1) is indeed closed to all orders in 
perturbation theory. Also there is no renormalization to A, 
whereas a correction to the elasticity J Zd x Z is allowed, and 
indeed shows up at 2-loop order. 

Finally, note that the domain of variation of u, in the peri- 
odic case yields an action with multiplicative periodicity in Z, 
but this does not seem to be important here. 

Let us now discuss the relation of our findings with self- 
avoiding polymers, branching processes and directed percola- 
tion. 

First, on figure IX. 4 we have drawn a diagram correspond- 
ing to the perturbation expansion of fixpoint SAP, which is 
the only fully attractive fixed point in the phase-diagram IX. 3. 
One easily checks that by integrating over times, one recovers 
a standard cf) 4 -perturbation theory, as depicted on figure IX.4. 
By first integrating over the momenta, one recovers the pertur- 
bation expansion of self-avoiding polymers. It is well known, 




FIG. IX.4: A self-avoiding polymer. 




a .'' 

time 



FIG. IX.5: Self-avoidance plus branching. 

that this fixed point is stable. In terms of particles, it can be 
interpreted as the world-lines of diffusing particles, which are 
not allowed to visit twice the same point in space. Let us now 
add some terms b and c. In interesting limit is A = oo, since 
there c can be set to zero. Adding a term propotional to b, the 
diffusing particle is allowed to branch. More precisely, two 
particles can meet at a time t. Then one of the particles be- 
comes inactive, before reappearing at some later time t' > t. 
One can interprete this as 

A + A — ► A + B (TX.17) 
B — ► A . (IX. 18) 

Particle B is completely inert, and does not diffuse away from 
its position of creation, before it decays into A again. How- 
ever note that any point in the future is equally likely to see 
B change back to A. This is very different from e.g. a spon- 
taneous decay. This process is depicted on figure IX.5. b can 
either come with a positive sign, or with a negative sign. If 
the sign is positive, this can be interpreted as the two particles 
attracting to make the branching-process. It is clear that after 
some critical threshold, the process and such the phase SAP 
becomes unstable, since the induced attraction between parti- 
cles tends to make them collapse at the same point in space 
and then annihilate. This leads to the runaway-flow in phase 
B-l on figure IX. 3. On the other hand, for negative b, even 
a large \b\ does not lead to a collapse. This is why on fig- 
ure VII. 1 in the case of A = oo the SAP-phase with a < 
extends to b — ► — oo. This remains valid for finite A if in 
the full flow-equations (VII. 6) to (VII. 8) c = is set from 
the beginning. However the situation for finite A discussed 
in (Vn.lO)-(VH.ll) maps in the language of branching pro- 
cesses to a finite initial ratio between a and b, parameterized 
by c = 6c~ A , which remains uncorrected under renormaliza- 
tion. The second branching-process c being present, it can 
render the phase SAP unstable to B-2. The vertex c is inter- 
preted as 

A — > C (IX. 19) 

A + C — ► A . (IX.20) 

This means that a particle A becomes spontaneously inactive 
at some time t. It remains at position x until at some time t' > 
t another particle A comes by to free it. The reduced flow- 
equations for the combined situation are given in (VII. 10) and 
(VII. 1 1), and lead to the instability of the phase SAP induced 
by the branching process c. 
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X. LONG-RANGE ELASTICITY 

Let us now study anisotropic depinning in the case of a man- 
ifold with long range (LR) elasticity, the elastic force in (II. 1) 
being, in Fourier: 

cq 2 u q j -» (c a \q\ a + cq 2 )u q , t (X.l) 

There are now two elastic constants, the LR one c a and the 
short range (SR) one c 2 , and we thus define the two dimen- 
sional regularization-parameters, 

e = 2a - d (X.2) 
K = 2-a. (X.3) 

The case of most interest corresponds to the parameters for 
the contact line depinning, d = 1, a = 1, i.e. e = k = 1. 

Power counting shows that disorder is perturbatively rele- 
vant below the critical dimension d < d c = 2a. Disorder is 
thus relevant for the contact line case but the crucial question 
we investigate here is whether the KPZ terms are important 
there. Study of the contact line depinning is usually performed 
within a d = 2a — e expansion (see Ref. [19]) at fixed a. This 
is the solid line in figure X. 1 . However as soon as elasticity is 
long range (k > 0) simple power counting shows that the KPZ 
terms are perturbatively irrelevant for d near d c . Working at 
fixed a as e.g. a = 1 is thus not the best method. One alter- 
native is to study the vicinity of the point d = 4, a = 2 and 
perform a double expansion both for e and n small. The idea is 
to determine a line c2kpz(<*) in the (a, d)-plane below which 
the KPZ terms are important and must be included. One can 
determine this line near the point d = 4, a = 2 and, by ex- 
trapolation, find on which side of the line lies the interesting 
case e = k = 1 (see figure X.l.). 

Through the replacement q 2 — > q a in the propagators of the 
1-loop diagrams of section V, it is easy to derive the 1-loop 
FRG equations for a general a, in presence of a KPZ term as 
in (II. 1). First one obtains as usual that c a is uncorrected to 
all orders, and thus we set c a = 1 in the following. Defining 
the dimensionless couplings 

A = AAJT C (X.4) 
A(«) - Af - e A(uAJ c ) (X.5) 

within a Wilson scheme where = Ac~ £ is the running UV 
cutoff, we find the flow equations: 

d e ln\ = (-K- A 2 A(0)- AA'(0+) (X.6) 
d e A(u) = (e - 2C)A(u) + uCA'(u) + 2X 2 A(u) 2 

-A'(u) 2 - A"(u)(A(u)- A(0)) . (X.7) 

We work to lowest order in both e and k (and thus neglect the 
small changes in the coefficients of order k) and define the 
ratios 

n = - , Ci = - ■ (X.8) 
e e 

Of course the SR part of the elasticity is corrected: 

8 e In (j^J = — k — AA'(0+) - A 2 A(0) (X.9) 



d 




1 2 

FIG. X.l: Phase-diagram in the (a, d)-plane. The solid line is e = 
0. The dashed line corresponding to (X.ll) separates the domain 
where an infinitesimal KPZ-term is relevant from those where it is 
irrelevant. At order e 2 this line will bend to the left, but should not 
cross the point d = 1 and a = 1. 

and we will focus on situations where it is irrelevant (a condi- 
tion which must be checked a posteriori). 

Note that since the LR-elasticity is uncorrected, the dimen- 
sionless variables, contrary to (V.5), are not divided by c 2 but 
by c a = 1 and their RG-equations thus do not contain addi- 
tional contributions from the corrections to c 2 . As a result A 
has now non-trivial corrections and the Cole-Hopf mapping 
no longer works, or has to be defined with a flowing A. 

Before embarking on a more detailed analysis let us indi- 
cate the main behaviour we expect from Eqs. (X.6) and (X.7). 
For A = one has the usual anisotropic depinning fixed point 
studied in Ref. [19]. One can perform a linear stability anal- 
ysis of this FP for small A. From (X.6) one finds that linear 
stability holds provided 

k > Ciso = | + 0(e 2 ) (X.10) 

for the non-periodic problem, and (i so = for the periodic 
case. This is the dashed line 

dKPz(a) = 5a-6 (X.ll) 

represented in Fig. X.l. For d > dKPz(a) the isotropic FP is 
stable. This is the case for the contact line depinning. On the 
other hand, one expects from Eqs. (X.6) and (X.7) that even 
then, if the value of A is large enough, the RG may flow again 
to KPZ strong coupling. This is the same run-away flow as 
for SR elasticity. Both fixed points should be separated by an 
instable fixed point, of which we will show that it is attainable 
perturbatively. Thus for d > dKPz(a) we expect, and find 
below, two phases: one where A flows to zero (denoted the 
ID phase) and one where the KPZ terms are important (the 
AD phase). The question is thus to determine the basin of 
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FIG. X.2: = (/t as a function of A for LR elasticity. Note that 
Ci(-A) = Ci(A). 



attraction of each phase and the critical (i.e. repulsive) fixed 
point which separates the two phases. Quite generally one 
expects a critical value A* below which A flows to zero and 
above which it runs away. 

A simple argument, confirmed by the more detailed analy- 
sis presented below, allows to estimate A* for small values of 
e, i.e. near d w cJkpzO*)' Since A (it) changes by order A 2 for 
small A, the Eq (X.6) gives the critical value: 



A* 



^ Ciso 

|A'(0+)| 



(X.12) 



for small k — (Ciso/ e )> where A'(0 + ) = O(e) takes its (nega- 
tive) value for the isotropic depinning fixed point. 

Although analysis of the full FRG-flow requires numerics, 
one can obtain some analytical information on the transition 
between the isotropic phase and the anisotropic strong-KPZ- 
coupling phase. 

A. Non-periodic systems 

Let us start with non-periodic systems and search for a pertur- 
bative fixed point of the system (X.6), (X.7). Interestingly in 
that case, there is one, whose properties depend continuously 
on k — k/c. 

For each value of k we can determine the FP through the 
following construction. Given the reparametrization invari- 
ance (V.6) of (X.7), we can always set 



A(0) = e 



(X.13) 



and for each fixed value of A search numerically for a fixed 
point function of Eq (X.7) which decreases at infinity (short 
range pinning force correlations of the random field type). 
Interestingly we find, through explicit numerical integration, 
that there is always one such solution, denoted by A! (it), 

if one tunes £ to a value noted C(A)- The resulting curve 
0(A) := C(A)/e is plotted in Fig. X.2. It starts at ((A = 
0) = 1/3 (the isotropic value) and increases as A increases. 



FIG. X.3: The function f(x) (green/light) and £i{x) (red/dark) de- 
fined in the text. One can read off £i as a function of k as follows: 
The curve y = f(x) yields A^ (x axis) from k (y axis) and in turn, 
one reads £i(k) (y axis) from A« (x axis) using the curve y — £i(x). 
This is indicated by the arrows. The result is plotted in Fig. X.4 



Considering the fixed point equation of (X.7) at u — using 
(X.13) shows that the value of A'(0 + ) is a simple expression: 



A!'(0 H 



eVl- 20(A) + 2A 2 . 



(X.14) 



Thus, reporting this value, as well as A(0) = e in Eq. (X.6) 
we see that for each value of k we can determine the value of 
A by solving the equation 



k = /(A) = Ci(A) - A + Ay 1 - 20(A) + 2A 2 . (X.15) 

Denoting A^ this solution, we obtain the FP function Ai (it) 

and the value of the roughness exponent O(^) := O(A^). 
Comparing (X.6) and (X.9) we note that the SR elastic part is 
indeed irrelevant as soon as £ > 0, and thus the above analysis 
is consistent. 

The curves /(A), and the resulting O(^) = 

C(Aft) = C(/ _1 («)) ^ plotted in Fig. X.3 and X.4 respec- 
tively. 




-1 -0.5 0.5 1 1.5 2 

FIG. X.4: Ci(k), using £(A) and A = 
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One sees that there is a solution with a positive X* k only if 
k > k c = -| consistent with the linear stability analysis given 
above. The roughness exponent associated to this FP then 
increases continuously, as shown on figure X.4, from Ci = | 
to larger values as k increases beyond k c . In particular, since 
we are interested in the point k = 1, e = 1 of the (a, d)-plane 
(see figure X.l) it is worth wile to give the extrapolation: 



C(k = 1) = 0.7 ; 



(X.16) 



and Aft=i = 1.037, values which give the simplest extrapola- 
tion for the contact-line depinning. One should however not 
expect too high a precision from this crude estimate. 

Thus we have found a non-trivial FP for this problem. It 
continuously depends on k/c and exists only for k/c > 1/3. 
The simplest scenario is that this FP is associated with the 
critical behaviour at the transition between the phase where 
KPZ is irrelevant (isotropic depinning) and the phase where 
KPZ grows (anisotropic depinning). To confirm it and check 
that this FP has only one unstable direction one needs a more 
detailed numerical analysis. Note that this is also indicated 
by an adiabatic approximation considering (X.6) alone and 
assuming that the disorder does not vary, which yields that the 
FP is repulsive if f'(x) > and attractive if f'(x) < 0. 

B. Periodic systems 

In the periodic case, since ( = is requested at any FP, we 
see that we cannot enforce the SR-elasticity coefficient c 2 to 
scale to under renormalization, since the FP condition on A 
implies that (X.9) vanishes. However if we start with a small 
ratio of c-i /c a or if the flow is such that this ratio gets small 
before we reach the fixed point, then it is legitimate to neglect 
the effect of C2. We restrict our analysis to that case, and study 
equations (X.6), (X.7) searching for a FP. A more detailed 
numerical analysis of the flow equation is feasable. 
It can easily be seen that the form 



A(u) = — (a + be 



- uX + ce uX ) 



(X.17) 



is not exactly preserved by the flow anymore (e.g. diA(u) 

yields a term proportional to ue~ uX di\i through variations 
of A; which here flows). One can still however search for 
exponential fixed points since then A; does not flow. (X.7) 
yields the conditions 



a + 2a 2 + Abe = 
b + Aab + b 2 + be = 
c + 4ac + c 2 +bc = 



(X.18) 
(X.19) 
(X.20) 



and we can set c = be A to ensure periodicity A(u) = A(l — 
u). We obtain the following fixed points: 



\f\ + 34c x + c 2 ~ x _ 
l + c x 



1 

"7 T 



4 Ay/l + Mc x + c 2X 



(X.21) 
(X.22) 



as well as two others (a = —\-,b = 0) and (a = 0, b = 0). 
The corresponding FP condition for A gives 



= — k — a — 2be 



(X.23) 



The FP with a positive b is the one of interest. It is again 
presumably the boundary between the zero and strong KPZ 
phases. The value of A^ is given by the positive root of 



4k = 1 + 



2 X -7 



\J\ + Me x + e 2 ~ x 



(X.24) 



which reproduces (X.12), A^ <~ 6k, to lowest order in k. 
One finds that A« increases monotonically with k and diverges 
Ak — ► +oo as k — > i. This suggests that for k > \ only the 
ID phase exist. 

XI. CONCLUSION 

In this paper we have reexamined the functional renormal- 
ization group approach to anisotropic depinning. This was 
mandatory since non-analytic renormalized disorder correla- 
tors were found to be crucial already for isotropic depinning 
and were neglected in previous approaches of AD. 

Indeed we have shown that the non-analyticity of disorder 
arising beyond the Larkin length is crucial to generate the 
KPZ-term, a first explicit field theoretic demonstration of how 
these terms appear at depinning. The resulting anomalous 
terms in the /3-function modify the flow compared to previous 
approaches in interesting ways. We found several non-trivial 
fixed points and for SR elasticity a Cole Hopf transformed the- 
ory which allows to simplify considerably perturbation theory 
and indicates that the KPZ coupling A/c is uncorrected to all 
orders. 

For LR-elasticity we have found the domains of parame- 
ters belonging to ID and AD respectively. We found that for 
the experimentally interesting case of contact-line depinning, 
two phases exist, ID and AD, and that the KPZ-coupling (i.e. 
the anisotropy) should be large enough for the AD class to 
apply (otherwise the ID exponents is expected [18, 26]). At 
the transition a larger value of ( w 0.7e (with e = 1 for the 
contact line) is obtained. This scenario could be checked in 
a numerical simulation. To make the comparison with exper- 
iments more accurate one should consider the more involved 
structure for the KPZ terms unveiled in [25] but this can be 
done by methods similar to the one introduced here. 

For SR-elasticity we have found interesting new fixed 
points. A bit disappointingly, they possess one unstable direc- 
tion and thus correspond to transient or critical behaviour, and 
not to the asymptotic behaviour which instead is controlled 
by a runaway flow to a regime not perturbatively accessible 
by the present method. On the other hand, an encouraging 
result is that we found a class of disorder correlators (in the 
form of exponentials) which should be invariant to all orders. 
These correspond to a set of branching processes which look 
tantalizingly close to the ones introduced to describe reaction 
diffusion and directed percolation. More work is necessary to 
understand this simpler equivalent class of theories at strong 
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coupling, as they may contain the key to this conjectured con- 
nection between anisotropic depinning and directed percola- 
tion (in d = 1 + 1) and its generalizations in terms of blocking 
surfaces (in higher d) and ultimately an understanding of the 
upper critical dimension for this problem. 

A posteriori, it is not suprising that the present approach 
yields again a flow to strong coupling KPZ, as it does in the 
thermal version of the problem [23, 24]. It is possible that 
as in the thermal problem another representation, as e.g. the 
directed polymer, better exposes the physics and in particular 
what is missed in the present approach. The corresponding 
formulation would be 



Z(x,t)= I V[y{ T )} 

ly{i)=x 

exp 



(XI. 1) 



i.e. a directed polymer in a random potential but with the 
choice T — \ and the additional self consistency condition: 



ny,r) = ^F^y,l]nZ(y,r) 



(XI.2) 



which relates the random potential to the pinning force and to 
the free energy of the directed polymer and makes the problem 
analytically far more complex. It may possess similar physics 
and thus be amenable to some extended FRG approach which 
would better account (as it does for the thermal problem) for 
the coarse grained correlations in the y direction a property 
clearly not taken into account by the present method, which 
treats correctly only correlations in the In Z space. 
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APPENDIX A: DIAGRAMS 

We use the following model setting c = rj = 1 to simplify 
notations 

S = rjiiu — cuAu — Au(Vu) (A.l) 

Jx.t 

-- u xt u xt > A(u xt - u xt >) (A.2) 

& Jx,t,t' 

One also has to specify a cut-off procedure. For convenience, 
we chose to put a mass-term. This is justified at 1-loop order 
since the results are universal, i.e. cut-off independent. At 
second order, one would have to be more careful and use, e.g. 
an external momentum IR cutoff. 

Many of the diagrams which we need are identical to the 
driven manifold problem at A = 0. These diagrams are de- 



tailed in [19]. The new diagrams are 



1 i 



= 2 [ f e ^( fc2 + m2 ) e - t '( fc2+m2 )|t-t'| 

J0<t<t' Jk 



x A'(0+)mt 
2XA'(0+)uu ( [ c -*'(fc 2 +™ 2 ) 



- ,,2, k 2 



(k 2 + m 2 ) 



(k 2 + m 2 ) 2 
{k 2 +m 2 )t') 



2AA'(0 + )ww 
k 2 



L 



k 2 k 2 
+ 



2(/c 2 + m 2 ) 3 {k 2 +m 2 y {k 2 +m 2 ) 3 
k 2 



AA' 



= AA'(0+)^. (fc2+ ^ 2)3 



>'(0+)A*/ 



(A.3) 




k 2 



{k 2 +m 2 f 



u(V«) 2 



(A.4) 



= 2A'(0+)A 



(k a +p a )p a 



k (0 + p) 2 + m 2 ){k 2 + m 2 ) 



'(0+)A f ■ 



P - 



2{kpf 



; 9A 

" ' (k 2 +m 2 ) 2 (k 2 + m 2 ) 3 

= r^ a . (0+)A / ; _^_„ a „ 

(A.5) 

(Note that A <-> — p 2 .) Dots indicate omitted subleading 
terms. 




-^A(0)A 3 
a .Ik 



(k 2 + m 2 ) 4 



(A.6) 



= -4A(0)A 2 / 
Jk 



k(k + p) (kp) 



-4A(0)A 2 

4 
d 



k (k 2 + m 2 ) 2 {(k+p) 2 +m 2 ) 



(kp) 2 2 _ k 2 (kp) 2 



-Xlwd 



k {k 2 + m 2 f (k 2 + m 2 ) 4 



m 2 ) 2 



(uAu) 
(A.7) 
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w(Vu) 2 (A.8) 




A(u xt - u xV ) A 



2 \2 I (k 2 ) 2 



I 



(A.9) 





(A. 10) 



Therefore, we have the following corrections to c, -q, A and 5 
(setting I := J k ( fc 2 4.^2)2 , dropping finite terms in e, but for 
the moment keeping the explicit d-dependence) 



5c/ c- 



2(2 rf) A'(0+)A-^A(0)A 2 
a a 



5r ] /? 1 = -\A'(0 + )I 

4aA'(0+) - 3A(0)A 2 + -U(0)A 2 



5X/X- 



d 



d 



d 



5A(u) = 2A{ u y\ 2 I 



(All) 
(A. 12) 
I 

"(A.13) 
(A. 14) 



The coupling-constant is A := A/c. Note that its flow vanishes 
at leading order in 1/e. We now check cancellations beyond 
the leading order. We use 



[ = A d r dkk^ 1 

Jk JO 

The two diagrams proportional to A'(0 + ) are: 



(A. 15) 



id 



k 2 



= A d 



k d{k 2 + m 2 ) 3 



(-2 + d) tt csc(^) 



4m e 



(A. 16) 




L 



Ak 2 



+ 



kd(k 2 + m 2 f (k 2 +m 2 ) 



-4, 



(2-d) tt csc(^) 
Am e 



(A47) 



The sum (which gives the renormalization of A) exactly van- 
ishes. 

The corrections proportional to A(0) are 



-8 k 




k d(k 2 +m 2 ) 4 



= A d 



-4+d 2 ) TT CSC(^) 

12 m e 



(A. 18) 




Ak 4 



k d(k 2 + m 2 Y 



= A d 



(A-d 2 ) TT CSC(^) 



24 m e 



(A. 19) 




4fc 2 



kd(k 2 +m 2 y d(k 2 +m 2 Y 



A (d-2 ) (d-l)TTCSc(^) 



12 m € 



The sum of the above three terms is 



A d 



Note that 



dn 

TT CSCl ) 

v 2 ; 



{A-d) (d — 2) 7r csc(^) 
24 m e 



7T 2 (d-4) 77T 4 (d-4) J 



(A.20) 



(A.21) 



d-4 ' 12 2880 + "' 

(A.22) 

So, working in a massive scheme, there are corrections at or- 
der e, compared to the leading term which would be 1/e. We 
see that the fixed point of Stepanow [22] is -even if one would 
accept his scheme- incorrect. However, as we have already 
stated above, one should do the calculations in a massless 
scheme. 

APPENDIX B: LONG RANGE DISORDER 

In this Appendix we give a quick study of the case with long 
range disorder in internal space x. We show that one recovers 
the Flory estimate of Section IV in the case of isotropic depin- 
ning. For anisotropic depinning we find a runaway flow and 
cannot conclude. 
We study 



Sbo = 77 



Z xt Z x > t >\ 2 : A(u x 



<~ x 



u x , v )f{x-x') (B.l) 
(B.2) 
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We find the FRG equation for the LR disorder: 

8A = eA + A(0)A" + (2- A1 )(AA'(0+) + A 2 A(0))A (B.3) 

with e = 4 — a, d large enough (d > a or more). We have 
absorbed eA in A with: 

A = f C(q) 2 f(q) (B.4) 
Jq 

This is because the graphs leading to two A (it) 2 functions or 
more do not contribute. This remains true to all orders, inspec- 
tion for A = shows that to two or three loops no corrections 
arise, except anomalous terms (which, as we will see are not 
needed as we find analytic fixed points). So for A = the one 
loop result is probably exact to all orders. 

The coefficient ^ comes from the corrections to the gradient 
term: 

(AA'(0+) + A 2 A(0))5 (B.5) 

B = hJ x2 f(^ c ^ (B - 6) 

fi = 2B/A= 2id d 4) (B.7) 

with a = 4, d > 4 (note that it goes to 2 when d goes to 
infinity). 

One easily finds the fixed points for A = 0. For periodic 
disorder one has: 

A(u) =.gcos(27rw) (B.8) 



dg = eg - (2n) 2 g 2 (B.9) 

The correlations are: 

( UU ) = A(0)^ ~ q-( d+2 V (B.10) 

with £ = e/2 as if A(0) was uncorrected. 
For non-periodic disorder, rescaling A gives: 

A(u) = A(0)e- e " 2 /( 6A (°» (B.ll) 
C = e/3 (B.12) 

Thus the Flory estimate is exact. Note that the fact that the LR 
correlator A(u) is analytic is not puzzling, since it generates 
in turn a SR part which should be nonanalytic in order to e.g. 
successfully generate a depinning threshold force. 

On the other hand for A > we find that: 

A(u) = ff cos(27rii) (B.13) 
dg = eg-(2TT) 2 g 2 + 1 X 2 g 2 (B.14) 

and there is thus a critical A beyond which there is no fixed 
point. This seems also to be the case for RF. Because of this 
runaway flow we cannot conclude. 
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